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Using a sample of 1.06 x 10* ip{2S) events collected with the BESIII detector at BEPCII, the 
decay 'ip{2S) — pprj is studied. A partial wave analysis determines that the intermediate state 
iV(1535) with a mass of 1524 ± Stf MeV/c^ and a width of HQtVit^ MeV/c^ is dominant in the 
decay; the product branching fraction is determined to be B{^{2S) — A''(1535)p) x _B(A'^(1535) 
prf) + C.C. = (5.2 ±0.3^^2) ^ 10~^. Furthermore, the branching fraction of ipi'^S) rjpp is measured 
to be (6.4 ±0.2 ±0.6) x 10"^ 
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I. INTRODUCTION 

Baryon spectroscopy is an important field to 
understand the internal structure of hadrons. Within 
the static quark model, the baryon octet and decuplct 
are well described. About half a century after 
the introduction of the quark model, however, a 
substantial number of light baryons predicted by the 
quark model have not been observed experimentally, 
which is known as the "missing baryon problem" 
0]. One possibility could be that the missing states 
simply do not exist, which has lead to the development 
of ncw_ phenomenological models, eg. the di-quark 
model [31]. Alternatively, the coupling of the unobserved 
states through conventional production channels could 
be small, which makes their observation more difficult. 

In addition to fixed target experiments [S-fTTj. 



charmonium decays produced in e+e" collisions open 
a window to hunt for the missing baryons [l^]- The 
Beijing Spectrometer (BES) [l^ experiment started a 
baryon program about a decade ago with the study of 
A^(1535) and iV(1650) in J/ip — >■ pp-q by partial wave 
analysis (PWA) using a sample of 7.8 million J/ip 
events. Using 58 million J/tp events collected at the 
BESII detector, a new excited nucleon 7V(2065) [H, 
was observed in J/^ pn-K~ (l7| and subsequently 
confirmed in J/ij: — >■ pp7r° [l^. BESII also studied 
'0(25') — )■ ppjj, where both ppir^ and ppi] were observed, 
ip{'2,S) — )■ pprj for the first time with a branching fraction 
of (5.8±1.1±0.7) X 10~^ . In both decays, there was weak 
evidence for a pp threshold mass enhancement but no 
PWA was performed [l^ . Most recently BESIII reported 
PWA results of ip{2S) — )■ ppn^ [l^l, and two new broad 
excited nucleons, 7V(2300) and Af(2570), were observed. 
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However, no clear evidence for 7V(2065) was found. 
Using 24.5 x 10^ ■0(25) events, CLEO-c [U reported 
the analysis of ip{2S) — jpp,Tr^pp and ripp without 
considering interference effects, in which A^(1535) and 
a pp enhancement(i?i(2100)) were investigated in ip{2S) 
decay to pprj. Those results show that J/ip and ip{2S) 
decays offer a unique place to study baryon spectroscopy. 

In this paper, using the 106 million ipi^S) events taken 
at the BESIII detector, a full PWA of the decay V'(2S') 
ppr] is performed. 



II. BESIII DETECTOR AND MONTE CARLO 
SIMULATION 



BEPCII [23 is a double-ring e+e collider designed 
to provide a peak luminosity of 10'^^ cm~^s~^ at a 
beam current of 0.93 A. The BESIII [Hi detector has 
a geometrical acceptance of 93% of in and consists of 
four main components: (1) A small-cell, helium-based 
(40% He, 60% CaHg) Main Drift Chamber (MDC) 
with 43 layers providing an average single-hit resolution 
of 135 /im, charged-particle momentum resolution in 
a 1 T magnetic field of 0.5% at 1 GeV/c^, and a 
dE/dx resolution, which is better than 6%. (2) A 
Time-Of-Flight system (TOP) constructed of 5-cm-thick 
plastic scintillators, with 176 detectors of 2.4 m length 
in two layers in the barrel and 96 fan-shaped detectors 
in the endcaps. The barrel (endcap) time resolution 
of 80 ps (110 ps) provides 2a K/tt separation for 
momenta up to ~ 1.0 GeV/c^. (3) An ElectroMagnetic 
Calorimeter (EMC) consisting of 6240 CsI(Tl) crystals 
in a cylindrical structure (barrel) and two endcaps. The 
energy resolution at 1.0 GeV is 2.5% (5%) in the barrel 
(endcaps), and the position resolution in the barrel 
(endcaps) is 6 mm (9 mm) . (4) The MUon Counter 
(MUC) consists of 1000 m^ of Resistive Plate Chambers 
(RPCs) in nine barrel and eight endcap layers and 
provides 2 cm position resolution. 

A GEANT4-based simulation software BOOST [H 
includes the geometric and material description of the 
BESIII detectors, the detector response and digitization 
models, as well as the tracking of the detector running 
conditions and performance. The production of the 
'0(25') resonance is simulated by the Monte Carlo 
(MC) event generator KKMC [2J], while the decays 
are generated by EvtCen [2^ for known decay modes 
with branching ratios being set to the PDG [2^ 
world average values, and by Lundcharm (27| for the 
remaining unknown decays. The analysis is performed 
in the framework of the BESIII Offline Software System 
(BOSS) [2^ which takes care of the detector calibration, 
event reconstruction and data storage. 



III. EVENT SELECTION 

For 0(25) ppr] (77 — > 77), the topology is quite 
simple, pfrf^. Each candidate event is required to have 
two good charged tracks reconstructed from the MDC 
with total charge zero. The point of closest approach 
to the bcamline of each charged track is required to be 
within ±20 cm in the beam direction and 2 cm in the 
plane perpendicular to the beam. Both tracks must have 
the polar angle 9 in the range of |cos0| < 0.93. The 
TOP and the specific energy loss dE/dx of a particle 
measured in the MDC are combined to calculate particle 
identification (PID) probabilities for pion, kaon and 
proton hypotheses. The particle type with the highest 
probability is assigned to each track. In this analysis, 
one charged track is required to be identified as a proton 
and the other as an anti-proton. 

Photon candidates are reconstructed by clustering 
EMC crystal energies. For each photon, the minimum 
energy is 25 MeV for barrel showers (| cos 6*1 < 0.80) and 
50 MeV for endcap showers (0.86 < jcos^l < 0.92). 
To exclude showers from charged particles, the angle 
between the nearest proton track and the shower must be 
greater than 10°, while for the anti-proton the angle has 
to be greater than 30°. Timing requirements are used 
to suppress electronic noise and energy deposits in the 
EMC unrelated to the event. At least two good photons 
are required. 

ii 1 
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FIG. 1: Scatter plots of pp invariant mass versus 77 invariant 
mass 

For the candidates remaining, a four-constraint 
kinematic fit imposing energy-momentum conservation 
is made under the ppjj hypothesis. If the number of 
selected photons is greater than two, the fit is repeated 
using all permutations of photons. The two photon 
combination with the minimum fit Xpp'yy is selected, and 
Xpp-yj is required to be less than 20. Because data and 
MC simulation do not agree well in the low momentum 
region, the momenta of the proton and anti-proton are 
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required to be greater than 300 MeV/c. Figure [T] 
shows the scatter plot of Mpp versus M^j for events 
satisfying the above requirements, where the two vertical 
bands correspond to the decays ip{2S) — )- ppn° and 
'0(25') — > ppr], and the horizontal band corresponds 
to the decay ■0(25') ^ X + J/0 (J/-0 pp). To 
remove the background events from V'(25) vJ/i' 
and 0^(25') -/Xcj, Mpp < 3.067 GeV/c^ and Mpp < 
(3.4 GeV/c^ — 0.75 x M^-y) are required. To select a 
clean sample, M^-y is required to be in the rj mass region, 
\M^^ - A/^l < 21 McV/c2. 

After the above event selection, 745 candidate events 
are selected. The Dalitz plot of Mpjj versus AIp^i is shown 
in Fig. [2] (a), where two clusters, corresponding to the pr] 
mass threshold enhancement displayed in Fig. [2] (b) and 
Fig. [5] (c) are visible. Both the mass spectra and the 
Dalitz plot display an asymmetry for prj and prj, which 
is mainly caused by different detection efficiencies for the 
proton and anti-proton. 

To investigate possible background events, the same 
analysis is performed on the MC sample of 100 million 
inclusive 0(25) events, and 11 backgound events are 
found from the channels, -0(25) — >■ JXcjiXcJ — >■ PP7i'°)i 
0'(25) -> -yxcjixcj iJ/i^, Jf-ip IPP) and 0(25) ^ 
7Xco(XcO — >■ pA+,A+ — > p7r°), which is compatible 
with the number of background events, 14, estimated 
with -q sidebands (|M^^ - 0.43| < 50 MeV/c^ and 
\M^^ - 0.65| < 50 MeV/c2). Additionally, 42.6 pb-^ 
of continuum data taken at 3.65 GeV/c^ is used for an 
estimation of the background from QED processes and 
51 background events are obtained after normalization 
with the luminosity of the continuum data and '0(25') 
data. The background events from 77 sidebands and 
the continuum data will be considered in the PWA of 
'0(25') — > pprj. 



IV. PARTIAL WAVE ANALYSIS 

The two-body decay amplitudes in the sequential decay 
process '0(25) — )■ N*p,N* — )• T]p (the charge-conjugate 
reaction is always implied unless explicitly mentioned) 
are constructed using the relativistic covariant tensor 
amplitude formalism [2^, and the maximum likelihood 
method is used in the PWA [H]. In 0(25) 
Nxp, Nx — > rjp, Aj is described as 



= Krod-xiBW)xAdecav-X] , 



(1) 



where -^prod-x ^^'^ amplitude, describing the 
production of the intermediate resonance Nx, BWx is 
the Brcit-Weigncr propagator of Nx, and Adecay-x is 
the decay amplitude of Nx- The total differential cross 



section ^ is 



phsp I 



(2) 



where Fphsp denotes the non-resonant contribution 
described by a interfering phase space term. The 
probability to observe the event characterized by the 
measurement ^ is 



P{0 = 



(3) 



where a;(^) = ^ and e(^) is the detection efficiency. 
/ d^uj{^)e{^) is the normalization integral calculated from 
the exclusive Monte Carlo sample. The joint probability 
density for observing n events in the data sample is 



/:=p(ei,6,-,u=n^(^»)=n 



i=l 



(4) 

Rather than maximizing the likelihood function ln(C), 
S = —InC is minimized to obtain Cj parameters, as well 
as the masses and widths of the resonances 



InC ■ 



EM 



c^(6) 



(5) 



i=l 



For a given data set, the second term is a constant and 
has no impact on the determination of the parameters 
of the amplitudes or on the relative changes of S values. 
So, for the fitting, —InC is defined as 



InC : 



EM 



Jd^ojiOeiO 



(6) 



The contribution of non-77 events and QED processes 
can be estimated with 77 sidebands and continuum 
data. In the log-likelihood calculation, the likelihood 
value of 77 sidebands and continuum data events are 
given negative weights, and are removed from data, 
since the log-likelihood value of data is the sum of the 
log-likelihood values of signal and background events 



5 = -[{-ln{C))data-{-ln{C)) 



bg\ 



(7) 



The free parameters are optimized by FUMILI [30| . In 
the minimization procedure, a change in log-likelihood of 
0.5 represents one standard deviation for each parameter. 

In the analysis, the following two Breit-Wigner 
formulas are used to describe the resonance. One form 
has a width which is independent of the energy of the 
intermediate state 



BW{s) 



1 



s — iMjv^Fjv. 



(8) 



where s is the invariant mass-squared. For Af(1535) with 
its mass close to the threshold of its dominant decay 
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FIG. 2: (a) The Dalitz plot of ip{2S) — >■ ?7pp and distributions for (b) M,,p, (c) Mr^p, and (d) Mpp. The crosses represent the 
data, the blank histograms show background events from continuum data, and the shaded histograms represent the background 
events from rj sidebands. 



channel Nt], the approximation of a constant width is 
not very good. Thus a phase space dependent width for 
Af(1535) is also used 



BW{s) 



1 



(9) 



s — iMis[*TM' (s) 
The phase space dependent widths can be written as [3l| 

r..(.). TV (0.5^^1^ + 0.5^^^), (10) 

where p^at and p^^r are the phase space factors for -kN 
and TjN final states, respectively, 



Pxn{s) 



2qxN{s) 

VI 

^{s - {Mn + MxY){s - {Mn - MxY) 



(11) 



V. SYSTEMATIC ERRORS 

The systematic error sources and their corresponding 
contributions to the measurement of mass, width and 
branching fractions are discussed below. 

• To investigate the impact on the PWA results 
from other possible components, the analysis is also 
performed including other possible N* states (eg. 
Af(1520), iV(1650), iV(1700), iV(1710), N{n2Q), 
iV(1895) and iV(1900)); the changes of the mass, 
width, and observed number of iV(1535) — > prj 
events are taken as the systematic errors by 
summing them in quadrature. 

• In the analysis, the background level is quite low, 
and the events from 77 sidebands and continuum 
data are considered in the PWA. To estimate 
the uncertainty, the background events from 77 
sidebands arew varied by ±50%, and the biggest 
change of the results is assigned as the systematic 
error. 



where X is tt or 77, and qxN{s) is the momentum of X in 
the center-of-mass system of XN . 



• In Eq. (|10p . the weight of the phase space factors 
for both ijN and ttN is set to be 0.5. The change 



IVl, p(GeV/c2) 



M JGeV/c^) 




Mp_(GeV/c2) cos(ep^) 

FIG. 3: Distributions of (a) Mp^, (b) Afp,,, (c) Mpp and (d) the angle between pr^ in the pp system. The crosses are for data, 
the blank histograms for PWA projections, the dashed lines for the contribution of A'^(1535) and the shaded histograms for the 
background events from -q sidebands and continuum data. 



of the results due to the variation of the weights in 
the range of '--^ 1 is taken as the systematic error. 

• The MDC tracking efficiency was studied with 
the clean sample of J/ip — >■ pptt~^tt~ events, as 
described in Rcf |32|. The difference between data 
and MC is less than 2% per charged track. Here, 
4% is taken as systematic error for the proton and 
anti-proton. 

• According to the particle identification efficiency 
study in Ref [s^l, the difference of the particle 
identification efficiencies between MC simulation 
and data is around 2% for each charged track. In 
this study, the two charged tracks are required to 
be identified as p and p, so 4% is taken as its 
systematic error from this source. 

• The systematic error from the photon detection 
efficiency has been studied using J/^ — > p°7r° 
events in Ref. [s^. The result indicates that 
the difference between data and MC simulation is 
about 1% for each photon. For the decay mode 
analyzed in this paper, 2% is taken as systematic 
error from two photons in the final states. 



In order to estimate the systematic error of the 
kinematic fit, a clean sample of J/ip — ppir'^ is 
selected. The difference of the efficiency between 
data and MC with and without using the four 
constraint kinematic fit, 7%, is taken as the 
systematic error. 



34|, 



• The number of ^/^ (25') events, (1.06±0.86) x 10® 
was determined from inclusive hadrons, and the 
systematic uncertainty is 0.82%. 

Table H] summarizes the systematic error contributions 
from different sources for the measurements of mass, 
width and branching fractions, and the total is the sum 
of them in quadrature. 



VI. RESULTS 

The PWA results, including the invariant mass spectra 
of pp, rjp, rjp and angular distributions, are shown as 
histograms in Fig. [3] and are consistent with the data. 
The best solution indicates that A^(1535) combined with 
an interfering phase space is sufficient to describe the 
data. We observe 527 ± 27 7V(1535) prj events 
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TABLE I: Summary of the systematic error contributions from different sources. 



Source 


A M(MeV/c2) 


A r(MeV/c2) 


A B(%) 


Additional Resonances 


+2 
-4 


+56 
-0 


+59 
-6 


Different BW Formula 






+17 
—21 


Background Uncertainty 




+10 
-10 


1 o 

+o 


MDC Tracking- 






±4 


Photon Detection 






±2 


Particle ID 






±4 


Kinematic Fit 






±7 


The number of '4>{'2S) events 






±0.82 


Total 


+10 

-4 


+57 
-10 


+63 
-23 



with a mass M = (1524 ± 51^'') MeV/c2, a width pprj is measured to be 



130 



+27+57 
-24-10 



MeV/c , and a statistical significance 
larger than lOcr. Here, the first error is statistical and 
the second is systematic. The contributions of Af(1535) 
and phase space are 70.8% and 61.0%. To determine the 
detection efficiency of ?/'(2S') A^(1535)p, the MC events 
are generated in accordance with the PWA amplitudes 
for V'(25') iV(1535)p. With the detection efficiency 
of 24.1%, the product branching fraction of ip{2S) — ?> 
7V(1535)p (7V(1535) pr]) is calculated to be 



B{^{2S) -> iV(1535)p (iV(1535) ^ pr/)) 

^obs n , n + 3.2 



£ ■ N^{2S) ■ B{t] -> 77) 



(5.2±0.3;j;2) X 10~^ 
(12) 



where the number of tp{2S) events, N^(2S)i is (1-06 ± 
0.86) X 10^ determined from '4>{2S) inclusive decays [s^ : 
B{-q ^ 77) is the world average value [111, and the first 
error is statistical and the second systematic. 

To investigate the pp mass enhancement observed at 
BESII [iil and CLEO-c [13, which did not use a PWA, 
a scan for an additional 1 resonance described by a 
Breit-Wigner function is performed. The widths used 
are 50 MeV/c^, 100 MeV/c^, 200 MeV/c^, 300 McV/c^, 
400 MeV/c2, 500 MeV/c^ and 600 MeV/c^. The mass 
is allowed to vary from 1900 MeV/c^ to 3000 MeV/c^ 
with steps of 2 MeV/c^. There is no evidence for a pp 
resonance in this region, indicating that the threshold 
enhancement can be explained by interference between 
the 7V(1535) and phase space. 

Subtracting the 51 and 15 background events from 
QED processes and from ry sidebands, respectively, the 
number of ^'(25') — > pprj events is calculated to be 
679 ± 26. In addition to the contribution from iV(1535), 
the contribution from the phase space events is taken into 
account in the determination of the detection efficiency 
according to the PWA results. With the detection 
efficiency of 25.6%, the branching fraction of ipi^S) — )• 



B{iP{2S) r)pp) = (6.4 ± 0.2 ± 0.6) X 10"^ (13) 



VII. SUMMARY 

Based on 1.06 x 10^ ■0(25') events collected with BESIII 
detector, a full PWA on the 745 ip{2S) — ppij candidates 
is performed, and the results indicate that the dominant 
contribution is from tp{2S) — > A^(1535)p. The mass 
and width of A^(1535) are determined to be 1524 ± 



5lf MeV/c2 and 130 



+27+57 
-24-10 



MeV/c , respectively, which 



are consistent with those from previous measurements 
listed in the PDG [35|. The product of the branching 
fractions is calculated to be B{ijj{2S) N{15S5)p) x 
B{N{1535) pr]) + c.c. = (5.2 ± 0.3t^^ J) x 10"^ The 
pp mass enhancement observed by BESII is investigated, 
and the statistical significance of an additional pp 
resonance is less than 3 <j. 

The branching fraction of 'il'{2S) — t- rjpp is determined 
to be (6.4±0.2±0.6) x 10^^, where the detection efficiency 
is determined from MC simulation events generated 
based on the PWA results. Compared with the branching 
fraction of J/i/j ppr] [11] , 



CPPV 



B{^j{2S) ^ Tjpp) 
B{J/iJ.' rjpp) 



(3.2 ±0.4)%, (14) 



which improves the BESII measurement [l0| of (2.8 ± 
0.7)%, and indicates that the decay ip{2S) — >■ pprj is 
suppressed compared with the "12% rule" . 
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